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In this work we calculate the decay rate of the electron type neutrinos into W bosons and electrons 
in presence of an external uniform magnetic field. The decay rate is calculated from the imaginary 
^ ' part of the W exchange neutrino self-energy diagram but in the weak field limit and compare our 

Q , result with the existing one. 
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I. INTRODUCTION 



The topic of neutrino propagation in various non-trivial medium has been widely studied. Particularly the topic 
'■^ ' dealing with the passage of neutrinos through the cosmos or through some compact object containing a magnetic field 
is of particular interest. As magnetic fields are ubiquitous in our universe so when neutrinos travel from their point 
O 1' of creation to their destiny they may feel the effect of the magnetic fields. As the neutrinos we consider are standard 
model chiral neutrinos they do not posses any magnetic moment and the magnetic field affects them through virtual 
, ■ charged particles inside the loops of the neutrino self-energy diagrams [3, 0, 9 • It is known that in these circumstances 
the magnitude of the magnetic field which is required to affect the neutrino properties is relatively big, O(10 13 )G. 
Such magnitude of magnetic fields are assumed to exist in neutron stars and magnetars. Consequently when neutrinos 
are produced in the neutron stars their properties can get modified because of the magnetic field. There has been 
many attempts to calculate the dispersion relation of neutrinos in presence of a magnetized medium, which employs 
the real part of the neutrino self-energy. But if we are trying to find out the absorption rate of the neutrinos then 
the imaginary part of the neutrino self-energy becomes important. There has been one attempt in the past [H where 
the authors tried to calculate the absorption coefficient of the neutrinos in a magnetic field background. In their 
work, Erdas and Lissia Q calculated the imaginary part of the self-energy of the neutrino, using the W exchange 
diagram, and calculated the absorption coefficient of the neutrino when it decays into an electron e~ and W~ pair, 
^— ^ in presence of a magnetic field. The present paper is aimed at calculating the damping rate of the neutrino due to its 
decay into e~ and W + pair, in presence of a magnetic field. The quantities named absorption coefficient by Erdas 
and Lissia and damping rate by the present authors turn out to be the same, both related with the attenuation of 
a traveling neutrino in presence of a magnetic field. Our work estimates the damping rate of neutrinos when the 
' neutrino propagation is damped due to its decay v e — > e~ + W + in presence of a magnetic field. The threshold of 
the pair production process v e — > v e + e~ + e + in presence of a magnetic field has also been calculated 0, @], and it 
is noted that the pair production threshold is smaller than the neutrino decay threshold. For higher energies of the 
neutrinos O(10 3_4 GeV) and for magnetic field magnitude O(10 n ~ 12 G) the neutrino decay rate becomes a prominent 
phenomena. The last reaction is kinematically forbidden in vacuum but can take place in presence of a magnetic field 
as in presence of a magnetic field the dispersion relation of the electron and positron changes. In their work Erdas 
and Lissia Q calculated the imaginary part of the neutrino self-energy using the Feynman gauge for the W boson 
propagator. In the present article we calculate the neutrino self-energy using the unitarity gauge. Our result is higher 
by 7% compared to the one by Erdas and Lissia. 

The paper is organized as follows: In sec. 2, we have shown the relation between the damping rate and the imaginary 
part of the neutrino self-energy. In Sec. 3 we explain the terminology of the Schwingcr proper time method and give 
the expression for the neutrino self energy due to W-cxchangc. The detail calculation of the imaginary part of the 
self-energy in the weak field limit and the damping rate are evaluated in Sec. 4. The threshold behavior as well as the 
asymptotic limit of the damping is also discussed in this section. A brief conclusion is given in Sec. 5. In appendix 
A we put some of the calculational details for the benefit of the reader. 
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II. DAMPING RATE OF CHIRAL FERMIONS 

In this section we briefly describe the main formula of the damping rate used in this paper. In general, the dispersion 
relation of a massless fermion with a momentum k^ = (k , k) is not given by k° = |k| and in particular, k is not zero 
at zero momentum. Writing 

k° = k° r -il, (2.1) 

where k® and 7 being real. The quantity M = k®(\k\ = 0) is interpreted as an effective mass and 7 as the fermion 
damping rate. With the above convention the probability of the chiral fermion not to decay in a time period t is 
~ e ft . The physically interesting regime is when 7 <C k®, since otherwise the system would be over damped and the 
concept of a propagating mode is not meaningful. If we decompose the self-energy of the chiral fermion in real and 
imaginary parts as 

S(fc) = ReE(fc) + ilmE(k) , (2.2) 

the damping rate comes out to be Q 

7 = --j-pU L ImE(k)u L , (2.3) 
where the spinors satisfy the Dirac equation 

(^-ReS(fc)K )fl = 0. (2.4) 

Eqn. (|2.3p can also be written as 

7 = -^oImTr[J(S](fc)], (2.5) 

where we treat the spinors as vacuum spinors and do neglect their corrections in the nontrivial media. In the present 
work we calculate the imaginary part of the self-energy of the standard model neutrino in presence of a strong magnetic 
field and using that result we find the damping rate of neutrino. To compare with earlier work we represent 7 in units 
of m _1 and not in the traditional s . 

III. THE W EXCHANGE DIAGRAM 

If the weak SU(2) coupling constant is g and the left-chiral projection operator is L = ^(1 — 75) then the contribution 
to the neutrino self-energy from the VK-boson exchange diagram is given by, 

-*S W r(fc) = J ^Li^^LiS e (p)(-i^=j ll/ LiW^(q) + c.t. (3.1) 

The c.t. is the counter term which has to be fed in the expression of the self-energy of the neutrino as the self-energy 
diagram contains ultraviolet divergence. The form of the c.t. is such that it cancels the divergence of the neutrino 
self-energy in the B — > limit. As the imaginary part of the self-energy in Eq. (|3.ip is relevant for the calculation of 
v e — > W + + e process so the c.t. will not affect our calculations as it does not contribute to the imaginary part of the 
self-energy. Consequently in this article we will not mention c.t. any further and write the self-energy of the neutrino 
in presence of a magnetic field as appearing in Eq. (|3.ip without the counter term. 

S e (p) is the electron propagator in presence of a magnetized medium and assuming the magnetic field to be along 
the 2-axis of the coordinate system S e (p)is given by: 

/>oo 

iS e {p)= / *»i e 0(p ' Sl) G(p,si). (3.2) 
Jo 

The possible phase term in the electron propagator is not written as it will disappear in the one-loop calculation of 
the neutrino self-energy. In the above equation, 

<f>(p,sx) = isx(po ~ ™ 2 ) ~ is i{pt+P±) tanZl ~ e|si I , (3.3) 

Zl 



W + (q = k-p) 




FIG. 1: The W exchange diagram. The solid internal virtual lines represent modified propagator in presence of an external magnetic 
field. 



where m is the mass of the electron and we have defined Z\ — |e|Bsi, p\ — p\ -\-p\ and e is an infinitesimal parameter 
introduced for the convergence of the integral in Eq. (|3.2p . Henceforth in this article we will not require the e term 
explicitly and so it will not appear in further discussions. For the sake of writing the electron propagator in a covariant 
form two 4- vectors, and b^, are used. As the magnetic field is a frame dependent quantity we assume in the rest 
frame of the observer u M is given by, 



(1,0) 



(3.4) 



We have a uniform magnetic field along the z-axis. Likewise the effect of the magnetic field enters through the 4- vector 
b 11 which is defined in such a way that the frame in which u M = (1,0) we have, 



V = (0, b) 



(3.5) 



where we denote the magnetic field vector by Bh. With the help of these two 4-vectors the other term appearing in 
Eq. (|3.2p can be written as, 



G(p, si) — sec 2 Z\ [A+ ipj5 + to(cos 2 z — iS 3 sinzi coszi)] 



(3.6) 



where and stand for, 



= Pp.- sin 2 2i (p ■ u Up- p ■ b 6 M ) , 
Bfj_ = sinzi coszi(p • u b^ - p ■ b u^) , 



(3.7) 
(3.8) 



and 



(3.9) 



The momentum dependence of the propagator are shown in Fig. [TJ The W boson propagator in presence of a magnetic 
field is W^ u {q) where q = k — p. In the weak-field limit i.e. eB <C M^, the VF-propagator is given by, 



(q 2 
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(3.10) 
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The q dependence of the terms in the denominator on the right hand side of the above equation can be expressed in 
integral form as 



1 



ds 2 e is ^~ M w) , 

) 

ds2S2 e - 2 (9 2 -M^) 



(3.11) 
(3.12) 



For convergence of the above integrals we require another term similar to the e term in Eq. (|3 - 2[) which has not 
been explicitly written as it will never appear in further discussions. Using Schwinger's proper time method the 
contribution to the neutrino self-energy coming from the W exchange process can be expressed as, 
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d A p 
d A p 



R^S e (p) lv LW^{q) 

pOO pOO 

i?7„ / d Sl / ds 2 e^^G(p, Sl ) lu L 



(2tt)4 

+ 2es 2 F» v ] e lS2(q2 - M ^ ) . 
As these integrals are lengthy and complicated to evaluate, we shall calculate them one by one. 



(3.13) 



IV. T, w (k) AND THE DAMPING RATE 7 

In this section we will evaluate T,w(k) given in Eq. (|3.13p . Before writing all the terms in Ew(fe) which will be 
evaluated separately at first we write down the phase appearing in the integrand on the right hand side of Eq. (13. 13[) . 
Let us define, 

A = <t>(p, si) + is 2 (q 2 - Ml,) 

z 2 fcii 2 z 1 k \ 2 

= -^( Sl m 2 + s 2 M^) + i(si + s 2 ){x 2 - x 2 ) - i\(x 2 + x\) - 1 j_ (4.1) 

zi + z 2 \e\B A \e\B 

where, 

kps 2 k 3 s 2 

xn=Po ; , 2:3 =P3 ■ , (4.2) 

si + s 2 si + s 2 

kis 2 k 2 s 2 
xi=pi — , x 2 =p 2 — , (4.3) 

Zi = \e\Bs l , A = 81 1- s 2 . (4.4) 

z\ 

Next the expression of T,w(k) m Eq. (|3.13[) is decomposed into various parts which are evaluated separately. The 
decomposition of S^(fc) is as follows: 



and 



IV(fc) = E#(fc) + E$(fc) + £$(fc) + E$(fc) , (4.5) 



where, 



■ 2 r j4 poc poo 

= \\ jAi J ds i J d S2 e A R 1 ,G(p,s 1 ) lv Lg^, 

■ 2 p j4 poo poo 



(4.6) 
(4.7) 
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4M 2 , J 


(2tt)4 
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(2^)4 



(4.9) 
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Doing the Gaussian integrals in the momenta we can write Eq. (|4.6p as, 



(4tt) 2 y 



(isi / cis 2 e A ° 



(si + s 2 )A 



i? 



S-2 



fc\\ + 



■S2 



(si + S 2 ) " («1 + S 2 ) 



-(kofi — k^fi) tan z\ 



•S-2 



where, 



sec zi 



A 



A: 



z 2 k± 



( Zl +z 2 )\e\B \\e\B \e\B 
The expression in Eq. (|4.7| can also be written as: 



(zim 2 + z 2 M 2 ) . 



(4.10) 



(4.11) 



■ 2 /* j4 poo poo 

S ^ (fc) = ~JM^Jj^ J, dSl J ^2e A sec 2 z 1 i?[2g.(A + zS) / 4-< Z 2 (A- 



= e£*>(*) + e£?>(*) 



where, 



(fc) = - 



'.9 



2M 2 , 



'.9 



d 4 p 
d 4 p 



2M 2 / J (2tt)4 y 



ds! / <fs 2 e A sec 2 z x .R [2g.(-4 + 1B)$ L . 
Jo 

/>0O 

dsi / ds 2 e A sec 2 zii? [<? 2 (A+ i$] L . 



(4.12) 

(4.13) 
(4.14) 



The elaborate expressions of the above integrals are given in Appendix [A"l The evaluation of S^/(fc) and S^/(fc) are 
straight forward and completing the Gaussian integrals over the loop momenta we obtain: 



2Af 2 (4tt) 2 J 



dsi I ds 2 e A ° 



S-2 



(ai + s 2 ) 2 A 



R 



i(ko J> — kj, jk) — tan z\ jku 



L. 



and 



2cBq f 00 f s ' 

Sy(fc) = . ,„ / dsi / ds 2 e A °- 2 R i(k fi- k 3 /i) -tan^i ^ii 

(4tt) 2 7 Jo (si + s 2 ) 2 A L 11 

The damping rate of the neutrino is given by, 



L . 



where 



-— Im«S(fe), 
2k 



S(k)=Tr[p: w (k)]. 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



In Eq. (|4.17[) we have used the free dispersion relation of the neutrinos to write the neutrino energy in the denominator 
and consequently here we have fcj? ~ k . The one-loop real part of the energy will differ from kg by factors proportional 
to the Fermi coupling Gf- For simplicity, we can change the integration variables from (si, s 2 ) to (z, u) where, 



m 1 

Z = (*1 + Z 2) = p(Zl + Z 2 ) , 



and 



M 2 , z 2 



1 z 2 



m z z\ + z 2 r)Zi + z 2 
With these new variables Eq. (|4.11|) can be written as: 



A 



o = -iz 



1 — UT\ + it + Wrj 



/If 2 

IV1 W 



1 



0Z 

7F 



(4.19) 



(4.20) 



(4.21) 
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where we have defined: 



Q = tan /3z(l — u) + (3zurj . 



(4.22) 



The phase exp(— iz) in exp(Ao) is oscillating and its main contribution to the z integral will come from z < 1. Also 
for weak-field limit the condition is eB <C m 2 . So in the weak field limit the product /?z< 1 and we can expand the 



terms in Yiy/(k) up to second order in (3z (Appendix [Al contains the expression of sj^- (fe)). This gives 



A — — iz 



1 — W] + It + 



Ik 2 , { eB 



u 2 z 2 (l - urif 



(4.23) 



(4.24) 



3 m 2 

As the £ l , for z = 1,2,3 and 4 are in the form: 

£ w (*) = /i%±) + /2(M-M). 

then from Eq. (|4. 18|) we see that the trace 

r &D£(*0)4-M)] = 0. (4.25) 
Consequently the terms which are proportional to Ln, will only contribute to the damping of the fcrmion. Also 

2 



noting that for a massless neutrino = and k 2 = k 2 ^ we can write S in the weak field limit, where /3z <C 1, as: 



S(k) 
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2 k 2 f°° f 1 ^ 



dz l du e 



4vr 2 Jo 

f3 2 r, 2 uz(u V - 1) {^2^4 _ ^02^ + u i(2z 2 f3 2 - + burj + 2} 



i/3 tj u{ut) — 1) 
90 



{z 2 j3 2 (ur] - l) 2 (100uV - 252«V + 129w?7 + 8) + 15(8uV - 22my + 11)} 



^From the above equation we can write, 



Here 



ImS(k) = ii^lH - 1 f dz [ duu {urj — 1) (T2 cos D + T\ sin D) 
12* J J 



Ti = z(2 + 5 ur) - u 2 r\ 2 ) + 2uV/3 * (! - 2mj + «V) , 

•Fa = ~r7(H - 22 u?7 + 8uV) + — ri{uT) - l) 2 /3 2 z 2 (8 + 129ur) - 252uV + 100uV) 
2 30 



(4.26) 



(4.27) 



and 



where we have defined 



By making the substitution 



we have 



D = (1 + u - wq)z + -£V(1 - urj) 3 z 3 , 



m \M? 



w 



' (£u) 2 / 3 (l-uri) 
D = xy H — y 3 , where x = 



1 q , (1 + u — urj) 



£2/3 M 2/ 3 (i - urj) 



(4.28) 



(4.29) 



(4.30) 



(4.31) 



(4.32) 
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In the new variable, ImiS(fc) can be written as: 



InuS 



12tt 2 ^S Jo 
In terms of the Airy function: 



dy 



1/77 



duu 3 



T 2 cos ( xy + -y 3 



sin j xy + -y 3 



1 f°° ( 1 . 

Ai(x) = - J dy cos \xy + -y' 



we can write Eq. (|4.33[) as: 



Im«S(fe) 



where 



12tt£§ 



du u 3 [(Wq - W 2 x + W 3 )Ai(x) - (Wi - W 3 x)Ai'(x)} 



W 
Wi 

w 2 
w 3 



-x, 



-77(11 - 22 urj + Su 2 ^ 2 ) . 

(2 + 5 urj - u 2 rj 2 ) 
1 + u — urj 

j](ur] - l) 2 /3 2 (8 + 129urj - 252«V + lOOuV ) 

30(l + u-tu7) 2 
2M 2 r; 2 /3 2 (l - 2 ur] + u 2 rj 2 ) 
(1 + it — W7) 3 



x 3 , 



(4.33) 
(4.34) 

(4.35) 

(4.36) 
(4.37) 
(4.38) 
(4.39) 



and A*'(aO = ^2. 



The numerical integration over u shows that, it saturates much before reaching the value 1 / 77, so here «ij < 1 and 
we can neglect the terms with urj or higher power of it in Eq. (|4.35p . So this gives 



ImS(fc) 



(jPvfkl f 1/v 

12<§ 



du u 3 



11?7 4 rj[3 2 



with 



! 15 (l+u) 2 ' 

(1 + u) 
(£u) 2 / 3 ' 



x A Ai(x) 



l + u 



Ai'(x) 



With these the damping rate of the neutrino is given as 
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V2G F k 2 ( eB 



M, 



H" 



1/77 



du u 3 
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/? 2 



1 2 15 Q • „: ••• 



l + u 



Ai'(x) 



(4.40) 



(4.41) 



(4.42) 



As 77 is much smaller than one we can henceforth neglect the part of the integrand proportional to 77 and assume 
I/77 ~ 00. Also here we assume that E v — ko ~ k± by taking k 3 to be very small and write the above equation as: 



7' 



V2. 
3tt 



--^—GpM w E u ( — 



2 / N 4 

TO 



r 2/3 / du ^_^ (a .) 
' 1+U 



(4.43) 



where S c = to 2 je is the critical magnetic field. For the transverse neutrino energy k± above the threshold of e and 
W production, we have observed that the function 



r 2/3 /(e) = ; 



-2/3 



du 



l + u 



Ai'{x) = 



-4/3 



du 



7^ J Jo 



2 (l + u) 3 / 2 



1. 



(4.44) 



In Eq. ()4.44j) , we express the derivative of the Airy function in terms of modified Bessel function to compare the u 
dependence of the integrand with Eq. (23) of ref.[4], which has a different u dependence and due to this our result 
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Er(eV) 



FIG. 2: The Eq. ()4.44p is plotted as a function of neutrino energy E v 
B/B c =0.1, 0.01 and 0.001 respectively. 



k± for different B. From left to right, the curves are for 



differs from the result by Erdas et al., in the above reference, 
damping rate behaves as 

i^g f m> w eJI 



For neutrino energy much above the threshold, the 



\Mw) 



(4.45) 



which is same as the absorption coefficient calculated by Erdas et al, in ref. [H . 

In Fig. 2, we have shown the behavior of Eq. (j4.44j) as a function of neutrino energy. For large k±, it saturates 
at 1.065 and below the energy threshold below the function £ _2 / 3 /(£) increases very rapidly and saturates very fast. 
Because of this behavior the damping rate also increases very rapidly from very small value and after crossing the 
threshold the behavior becomes linear in neutrino energy. For small magnetic field, the saturation energy is larger 
than the one with larger magnetic field. As can be seen from the figure, for B/B c =0.1, the curve saturates around 
k± ~ 10 19 eV, whereas for B/B c =0.001 the energy is about 10 21 eV. For the neutrino energy much above the threshold 
the asymptotic behavior is given in Eq. (|4.45[) . where the damping rate is proportional to the neutrino energy and 
also it depends quadratically on the magnetic field as shown in the above equation. The behavior can be seen from 
Fig. 3, where we have plotted the damping rate as a function of energy for three different values (B/B c = 0.1, 0.01 
and 0.001) of the magnetic field. Again larger the magnetic field, larger is the damping rate and for neutrino energy 
smaller than the threshold one, the damping rate is suppressed but increases rapidly by increasing the energy. 

The behavior of the process v e — ► e~ + W + solely depends on £~ 2 / 3 /(i;) and the threshold condition satisfies 
£~ 2 ^ 3 /(£) < 1- We have analyzed the threshold behavior of the above process by taking into account the behavior of 
the functions £~ 2 / 3 , /(£) and their product as function of neutrino energy, which are shown in Fig. 4 for B/B c = 0.1 
and accurately find out the threshold energy. With increasing energy the function £~ 2 / 3 decreases and / (£) increases 
and both of them intersect at a point ~ 0.66 which lies before the saturation point of their product and this is the 
point of threshold. The estimated neutrino energy at this point is E v — 2.37 x 10 17 eV. We have observed that, for 
B/B c = 0.01 and 0.001, the intersection point of the two functions as described above remains the same (0.66) but 
the neutrino energy shift towards higher values. For B/B c = 0.01, we obtain E v — 2.37 x 10 18 eV and 2.37 x 10 19 eV 
respectively. For v„ and v T decay calculations the electron mass has to be replaced by the corresponding lepton mass. 



1 = I>7 '.E 



(4.46) 
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FIG. 3: The damping rate is expressed in units of meter inverse is plotted against the transverse energy of the neutrino E y . From 
the top to bottom are the curves for B/B c = 0.1, 0.01 and 0.001. 




FIG. 4: The functions £ 2/s , /(£) and £ 2/3 /(£) are plotted as function of E v for J3/B c = 0.1. The extreme left decreasing curve 
is the function £ _2//3 , the second increasing curve is /(£) and the third one which is saturated around 1 is the function £~ 2 / 3 /(£). 
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or in terms of the above, we can write the threshold neutrino energy as 

Ef = 2.36 x 10 16 (j^J eV (4.47) 
which is about 7% higher than the one obtained by Erdas et. al.j4|. 

V. CONCLUSION 

By using the Schwinger's proper time electron propagator we calculate the imaginary part of the neutrino self 
energy in a constant magnetic field background. In this calculation we assume that the magnetic field to be weak 
eB -C m 2 . From the imaginary part of the neutrino self energy we calculate the damping rate i.e. conversion of 
neutrino into electron and W-boson. We explicitly evaluated all the contributions in the weak field limit. We found 
that the behavior of the process v e — > e~ + W + solely depends on the quantity £~ 2 / 3 /(£)- The threshold energy 
for this process can be accurately determined from the intersection of the functions £~ 2 / 3 and /(£) as functions of 
energy. Also for the neutrino energy much above the threshold we found that the damping rate is proportional to 
the neutrino energy and also it depends quadratically on the strength of the magnetic field. Our calculation gives an 
alternative way to arrive at the absorption cross-section calculated previously in 0]. The damping rate expressed in 
dimensions of inverse length matches to a reasonable degree with the old result although the gauge to fix the weak 
boson propagator is different in our case. This shows the gauge invariance of the calculations also. 



APPENDIX 



APPENDIX A: THE INTEGRAL REPRESENTATION OF THE W EXCHANGE DIAGRAM 



The expression in Eq. (|4.7[) can also be written as: 

.2 r rfi poo /-oo 

J (2^4 J dSl J ^ s 2 e A sec 2 Z\R [2q.(A + iB) fa - q 2 (A+ ^)] L , 



19 
2M 2 , 



where, 



HI 



2M 2 , 



dsi / ds 2 e. sec 2 z x R [2q.(A + iB)fa] L , 



2M 2 , 



J d Sl J ds 2 e A sec 2 z x R [g 2 (A+ 1% L . 



(Al) 

(A2) 
(A3) 



The term q.(A + iB)fa in the integrand on the right hand side of Eq. (|A2[) can be written in terms of xq, x\, x 2 , £3 as, 



q.(A + iB)i = 



d — XqCOS Z\ + XqXq + X\ + x\ + x\ COS Z\ + 



(xifci + x 2 k 2 ) 
A 



s 2 - Sl- 



tan Z\ 



- A3X3 



x (oq - 7o£o +7i xi + "f 2 x 2 + 73X3) . 



where, 



and, 



a 



si + S2 



si tan z\ 
A z\ 



d 



sis 2 2 2 sis 2 tanzi 2 

rfell COS Z\ K± . 



(*i + s 2 ) 2 



A 2 



An = 



(si - s 2 ) 
(si + s 2 ) 



ko cos 2 Z\ + i sin z\ cos Z]k? >1 A3 = 



(si - s 2 ) 
(si + 82) 



k$ cos 2 Z\ + i sin z\ cos Zi fen. 



(A4) 
(A5) 

(A6) 



where we have assumed that fcn ~ k± by taking k% to be very small. Expressing the derivative of the Airy function 
in terms of Bcsscl function as 



(A7) 
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we can express the damping rate as 



y/2 1 
V2 



'„ 2 /3 



G, 



-(\e\B) 2 ^(m 2 M w ) 2 ^ 



du 



173^2/3 ^3 — 



2 (1 + u) 3 / 2 



= ^-G F M 2 k ± — 



3tt 



D \ 2 / \ 4 

B \ / TON 



Bj \M W J Vv^ 4 / 3 



1/^2/3 — 



2(l + u)3/ 2 



In a similar fashion the term q 2 (A+ in Eq. (|A3j) can be written as, 



? 2 (A+^) 



I Z A Z A 

yo T x ~ *1 ~ x 2 ~ x 3 



(si + s 2 ) 

X (h + h Q XQ - 71X1 - 72X2 - h 3 x 3 ) , 



2si 2si tanzi 

a^ofco H — ; (a;ifci + X2K2) 



2si 



where, 



,9o 

h 



2s 2 



(si + s 2 ) 2 (si + s 2 ) 



h 2 ' 



A z\ 



2s 2 



Sl + S 2 



A 2 A 



jk\\ cos 2 z\ ^- k, + i — sinzi coszi(ko¥> — k^fy 

A S\+ S2 



■si + s 2 



and 



/lo = cos 2 z\t/x+ zJ6sinzi coszi, ft. 3 — cos 2 z\^> + ijlsinzi coszi 



(A8) 



(A9) 

(A10) 
(All) 

(A12) 



Using the expressions of q.{A + iB)<f[ and q 2 (A+ ifi) in Eqns. (|A2|) and (|A3|) and doing the integrals we obtain: 



^(2a) 



(*) = 



M 2 , (4^) 2 7 



dsi / <is 2 e 



(si + s 2 )A 



x R 



S-i S2 9 2 

fcii COS Zi 



SfS 2 



tan z\ \ 



+i 
9 2 



ll U(*i + s 2 ) 3 1 

/ s 2 S2 tan 2 Z\ 



(si+s 2 )A 2 x zi 
s 2 s 2 sin^! coszi 



2si 



(3si - s 2 ) 2 

cos zi 



(si + s 2 )A (si + s 2 ) 



(si +s 2 ) 2 A 



2si sinzi coszi 1 tanzi . \ ") sinzi coszi 

- + T^(3si S 2 ) )> + (fc 3 )4- fc p) ; 

A 2 z\ I I si + s 2 



(si + s 2 )A zi 
1 



(A13) 



2M 2 , (4tt) 2 y 



dsi / <is 2 e 



(si + s 2 )A 



x R 



h 



(2si + s 2 ) , 2 



'IV (^i + s 2 ) 3 11 (si + s 2 )A 2 



tan Z\ 



2si ' fci 

21 



+(*o]A-M 



tan zi 2 

K 1 



21 



s 2 



Sl + S 2 



■ sm z\ cos z\ 



s 2 (2ai +s 2 ) 
A (si + s 2 ) 2 
2si 



cos 2 zi — 2i cos Zi 
2s 2 



si - s 2 



s 2 



(si + s 2 ) 2 (si+s 2 )A 
2s 2 2sitanz!\ 



(si + s 2 )A A 2 A zi 
2 2^ 



s 2 (2si + s 2 ) 2 2sis 2 tanzi 2 

fZ\\ " — ~ ^ K I 



(si + s 2 ) 



A 2 



S2(si+s 2 ) (si + s 2 ) A 



-21 



(A14) 



(2) 

This results are used to evaluate T, w (k) 
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